CE-QUAL-W2 Data Visualization Demo

import w2

Meteorology: Load, diplay, and plot the data

data_columns = [
    'Air Temperature',
    'Dew Point Temperature',
    'Wind Speed',
    'Wind Direction',
    'Cloudiness',
    'Solar Radiation'
]

ylabels = [
    'Air Temperature ($^oC$)',
    'Dew Point Temperature ($^oC$)',
    'Wind Speed (m/s)',
    'Wind Direction (radians)',
    'Cloudiness (fraction)',
    'Solar Radiation ($W/m^2$)'
]

met_wb1 = w2.read('BerlinMilton2006/2006_Met.npt', 2006, data_columns)
met_wb2 = w2.read('BerlinMilton2006/2006_Met1.npt', 2006, data_columns)

display(met_wb1)
display(met_wb2)

w2.multiplot(met_wb1, ylabels=ylabels, colors=w2.k2)
w2.multiplot(met_wb2, ylabels=ylabels, palette='colorblind')

	Air Temperature	Dew Point Temperature	Wind Speed	Wind Direction	Cloudiness	Solar Radiation
Date
2006-01-01 00:00:00	0.6	-0.6	3.89	4.01	10.0	0
2006-01-01 01:00:00	0.6	-0.6	3.89	4.01	10.0	0
2006-01-01 02:00:00	1.1	-0.6	3.26	4.01	10.0	0
2006-01-01 03:00:00	1.1	-0.6	3.96	4.04	10.0	0
2006-01-01 04:00:00	1.1	-0.6	3.58	4.16	10.0	0
...	...	...	...	...	...	...
2006-12-31 20:00:00	11.0	7.4	7.09	2.44	10.0	0
2006-12-31 21:00:00	10.6	8.5	7.04	2.44	10.0	0
2006-12-31 22:00:00	11.1	9.5	8.34	2.47	10.0	0
2006-12-31 23:00:00	11.3	9.6	5.12	2.80	10.0	0
2007-01-01 00:00:00	11.3	9.6	5.12	2.80	10.0	0

8761 rows × 6 columns

	Air Temperature	Dew Point Temperature	Wind Speed	Wind Direction	Cloudiness	Solar Radiation
Date
2006-01-01 00:00:00	0.6	-0.6	3.89	4.01	10.0	0
2006-01-01 01:00:00	0.6	-0.6	3.89	4.01	10.0	0
2006-01-01 02:00:00	1.1	-0.6	3.26	4.01	10.0	0
2006-01-01 03:00:00	1.1	-0.6	3.96	4.04	10.0	0
2006-01-01 04:00:00	1.1	-0.6	3.58	4.16	10.0	0
...	...	...	...	...	...	...
2006-12-31 20:00:00	11.0	7.4	7.09	2.44	10.0	0
2006-12-31 21:00:00	10.6	8.5	7.04	2.44	10.0	0
2006-12-31 22:00:00	11.1	9.5	8.34	2.47	10.0	0
2006-12-31 23:00:00	11.3	9.6	5.12	2.80	10.0	0
2007-01-01 00:00:00	11.3	9.6	5.12	2.80	10.0	0

8761 rows × 6 columns

QIN (inflows), Berlin Model, 2006

Input files:

• BR1: 2006_DeerCrk_Qin.npt
• BR2: 2006_WillowCrk_Qin.npt
• BR3: 2006_IslandCrk_Qin.npt
• BR4: 2006_MillCrk_Qin.npt

data_columns = ['Inflow (QIN)']

qin_br1 = w2.read('BerlinMilton2006/2006_DeerCrk_Qin.npt', 2006, data_columns)
qin_br2 = w2.read('BerlinMilton2006/2006_WillowCrk_Qin.npt', 2006, data_columns)
qin_br3 = w2.read('BerlinMilton2006/2006_IslandCrk_Qin.npt', 2006, data_columns)
qin_br4 = w2.read('BerlinMilton2006/2006_MillCrk_Qin.npt', 2006, data_columns)

ylabel = 'Inflow ($m^3/s$)'
w2.plot(qin_br1, ylabel=ylabel)
w2.plot(qin_br2, ylabel=ylabel)
w2.plot(qin_br3, ylabel=ylabel)
w2.plot(qin_br4, ylabel=ylabel)
qin_br1

	Inflow (QIN)
Date
2006-01-01 00:00:00	1.02
2006-01-01 01:00:00	0.92
2006-01-01 02:00:00	0.91
2006-01-01 03:00:00	1.00
2006-01-01 04:00:00	1.00
...	...
2006-12-31 20:00:00	1.07
2006-12-31 21:00:00	0.98
2006-12-31 22:00:00	1.07
2006-12-31 23:00:00	1.26
2007-01-01 00:00:00	1.36

8761 rows × 1 columns

QOT (dam outflows), Berlin Model, 2006

Two outlets at Berlin Dam

Input files:

• BR1: 2006_BerlinDam.npt
• BR2: 2006_MiltonDam.npt

data_columns = ['Dam Outflow (QOT): Outlet 1', 'Dam Outflow (QOT): Outlet 2']

qot_br1 = w2.read('BerlinMilton2006/2006_BerlinDam.npt', 2006, data_columns)
qot_br2 = w2.read('BerlinMilton2006/2006_MiltonDam.npt', 2006, data_columns)

ylabel = 'Dam Outflow ($m^3/s$)'
w2.plot(qot_br1, title='2006_BerlinDam.npt', ylabel=ylabel)
w2.plot(qot_br2, title='2006_MiltonDam.npt', ylabel=ylabel)

QTR (tributary inflows), Berlin Model, 2006

Input files:

• TR1: 2006_BeechCrk_Qtr.npt
• TR2: 2006_Alliance_Qtr.npt
• TR3: 2006_MahoningR_Qtr.npt

data_columns = ['Tributary Inflow (QTR)']

qtr_tr1 = w2.read('BerlinMilton2006/2006_BeechCrk_Qtr.npt', 2006, data_columns)
qtr_tr2 = w2.read('BerlinMilton2006/2006_Alliance_Qtr.npt', 2006, data_columns)
qtr_tr3 = w2.read('BerlinMilton2006/2006_MahoningR_Qtr.npt', 2006, data_columns)

ylabel = 'Tributary Inflow ($m^3/s$)'
w2.plot(qtr_tr1, title='2006_BeechCrk_Qtr.npt', ylabel=ylabel)
w2.plot(qtr_tr2, title='2006_Alliance_Qtr.npt', ylabel=ylabel)
w2.plot(qtr_tr3, title='2006_MahoningR_Qtr.npt', ylabel=ylabel)

TTR (tributary inflow temperature), Berlin Model, 2006

Input files:

• TR1: 2006_MahoningR_Ttr5.npt
• TR2: 2006_Alliance_Ttr.npt
• TR3: 2006_MahoningR_Ttr.npt

data_columns = ['Inflow Temperature (TTR)']

ttr_tr1 = w2.read('BerlinMilton2006/2006_MahoningR_Ttr5.npt', 2006, data_columns)
ttr_tr2 = w2.read('BerlinMilton2006/2006_Alliance_Ttr.npt', 2006, data_columns)
ttr_tr3 = w2.read('BerlinMilton2006/2006_MahoningR_Ttr.npt', 2006, data_columns)

ylabel = 'Inflow Temperature ($^oC$)'
w2.plot(ttr_tr1, title='2006_MahoningR_Ttr5.npt', ylabel=ylabel)
w2.plot(ttr_tr2, title='2006_Alliance_Ttr.npt', ylabel=ylabel)
w2.plot(ttr_tr3, title='2006_MahoningR_Ttr.npt', ylabel=ylabel)

CTR (tributary constituent inflow concentrations), Berlin Model, 2006

Input files:

• TR1: 2006_MahoningR_Ctr1.npt
• TR2: 2006_AllianceWwtp_Ctr.npt
• TR3: 2006_MahoningR_Ctr.npt

infiles = [
    '2006_MahoningR_Ctr1.npt',
    '2006_AllianceWwtp_Ctr.npt',
    '2006_MahoningR_Ctr.npt'
]

ctr1_vars = ['TDS', 'SO4', 'Cl', 'ISS', 'OP', 'NH4', 'NOx', 'Fe',
             'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']
ctr2_vars = ['TDS', 'ISS', 'OP', 'NH3', 'NOx', 'LDOM', 'RDOM', 'LPOM', 'RPOM', 'DO']
ctr3_vars = ['TDS', 'SO4', 'Cl', 'ISS', 'OP', 'NH4', 'NOx', 'Fe',
             'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']

ctr1_ylabels = [
    'TDS (mg/L)', 'SO4 (mg/L)', 'Cl (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH4 (mg/L)', 'NOx (mg/L)', 'Fe (mg/L)',
    'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']

ctr2_ylabels = ['TDS (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH3 (mg/L)', 'NOx (mg/L)',
                'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'DO (mg/L)']
ctr3_ylabels = [
    'TDS (mg/L)', 'SO4 (mg/L)', 'Cl (mg/L)', 'ISS (mg/L)', 'OP (mg/L)', 'NH4 (mg/L)', 'NOx (mg/L)', 'Fe (mg/L)',
    'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']

ctr1 = w2.read(f'BerlinMilton2006/2006_MahoningR_Ctr1.npt', 2006, ctr1_vars)
ctr2 = w2.read(f'BerlinMilton2006/2006_AllianceWwtp_Ctr.npt', 2006, ctr2_vars)
ctr3 = w2.read(f'BerlinMilton2006/2006_MahoningR_Ctr.npt', 2006, ctr3_vars)

w2.multiplot(ctr1, ylabels=ctr1_ylabels, figsize=(15, 30), style='.-')
w2.multiplot(ctr2, ylabels=ctr2_ylabels, figsize=(15, 19), style='.-')
w2.multiplot(ctr3, ylabels=ctr3_ylabels, figsize=(15, 30), style='.-')

QDT (distributed tributary inflows), Berlin Model, 2006

Input files:

• BR1: 2006_Bqwb_72smooth.npt
• BR5: 2006_Mqwb_72smooth.npt

data_columns = ['Distributed Tributary Inflows (QDT)']

qdt_br1 = w2.read('BerlinMilton2006/2006_Bqwb_72smooth.npt', 2006, data_columns)
qdt_br5 = w2.read('BerlinMilton2006/2006_Mqwb_72smooth.npt', 2006, data_columns)

w2.plot(qdt_br1, title='2006_Bqwb_72smooth.npt', xlabel='Date', ylabel=f'Distributed Tributary Inflows ($m^3/s$)')
w2.plot(qdt_br5, title='2006_Mqwb_72smooth.npt', xlabel='Date', ylabel=f'Distributed Tributary Inflows ($m^3/s$)')

TDT (distributed tributary inflow temperatures), Berlin Model, 2006

Input files:

• BR1: 2006_Btdt.npt
• BR5: 2006_Mtdt.npt

data_columns = ['Distributed Tributary Inflow Temperature (TDT)']

tdt_br1 = w2.read('BerlinMilton2006/2006_Btdt.npt', 2006, data_columns)
tdt_br5 = w2.read('BerlinMilton2006/2006_Mtdt.npt', 2006, data_columns)

w2.plot(tdt_br1, title='2006_Btdt.npt', ylabel=f'Distributed Tributary Inflow Temperature ($^oC$)')
w2.plot(tdt_br5, title='2006_Mtdt.npt', ylabel=f'Distributed Tributary Inflow Temperature ($^oC$)')

CDT (distributed tributary constituent inflows), Berlin Model, 2006

Input files:

• BR1: 2006_MahoningR_Ctr1.npt

• BR5: 2006_MahoningR_Ctr2.npt

  Note: We don't need to analyze these because the CDT inputs use the CTR input files, which have been analyzed above

PRE (precipitation), Berlin Model, 2006

Input files:

• BR1: 2006_Pre.npt
• BR2: 2006_Pre1.npt
• BR3: 2006_Pre2.npt
• BR4: 2006_Pre3.npt
• BR5: 2006_Pre4.npt

data_columns = ['Precipitation']

pre_br1 = w2.read('BerlinMilton2006/2006_Pre.npt', 2006, data_columns)
pre_br2 = w2.read('BerlinMilton2006/2006_Pre1.npt', 2006, data_columns)
pre_br3 = w2.read('BerlinMilton2006/2006_Pre2.npt', 2006, data_columns)
pre_br4 = w2.read('BerlinMilton2006/2006_Pre3.npt', 2006, data_columns)
pre_br5 = w2.read('BerlinMilton2006/2006_Pre4.npt', 2006, data_columns)

ylabel = 'Precipitation (mm)'
w2.plot(pre_br1, title='2006_Pre.npt', ylabel=ylabel)
w2.plot(pre_br2, title='2006_Pre1.npt', ylabel=ylabel)
w2.plot(pre_br3, title='2006_Pre2.npt', ylabel=ylabel)
w2.plot(pre_br4, title='2006_Pre3.npt', ylabel=ylabel)
w2.plot(pre_br5, title='2006_Pre4.npt', ylabel=ylabel)

TPR (precipitation temperature), Berlin Model, 2006

• BR1: 2006_PreT_noneg.npt
• BR2: 2006_PreT_noneg1.npt
• BR3: 2006_PreT_noneg2.npt
• BR4: 2006_PreT_noneg3.npt
• BR5: 2006_PreT_noneg4.npt

data_columns = ['Precipitation Temperature (TPR)']

tpr_br1 = w2.read('BerlinMilton2006/2006_PreT_noneg.npt', 2006, data_columns)
tpr_br2 = w2.read('BerlinMilton2006/2006_PreT_noneg1.npt', 2006, data_columns)
tpr_br3 = w2.read('BerlinMilton2006/2006_PreT_noneg2.npt', 2006, data_columns)
tpr_br4 = w2.read('BerlinMilton2006/2006_PreT_noneg3.npt', 2006, data_columns)
tpr_br5 = w2.read('BerlinMilton2006/2006_PreT_noneg4.npt', 2006, data_columns)

ylabel = 'Precipitation Temperature ($^oC$)'
w2.plot(tpr_br1, title='2006_PreT_noneg.npt', ylabel=ylabel)
w2.plot(tpr_br2, title='2006_PreT_noneg1.npt', ylabel=ylabel)
w2.plot(tpr_br3, title='2006_PreT_noneg2.npt', ylabel=ylabel)
w2.plot(tpr_br4, title='2006_PreT_noneg3.npt', ylabel=ylabel)
w2.plot(tpr_br5, title='2006_PreT_noneg4.npt', ylabel=ylabel)

---

Analyze Outputs

---

CWO (constituent outflows), Berlin Model, 2006

• SEG 37: cwo_37_wdo.csv
• SEG 75: cwo_75_wdo.csv

# Constituents list: Include an extra "constituent" to the extra column that will be added due to a trailing comma.
constituents = ['TDS', 'G1_SO4', 'G2_Cl', 'ISS1', 'PO4', 'NH4', 'NO3', 'FE', 'LDOM', 'RDOM', 'LPOM', 'RPOM', 'BG', 'DIAT', 'OTH', 'DO']
cwo_37 = w2.read('BerlinMilton2006/cwo_37_wdo.csv', 2006, constituents)
cwo_75 = w2.read('BerlinMilton2006/cwo_75_wdo.csv', 2006, constituents)

ylabels = [
    'TDS (mg/L)', 'G1_SO4 (mg/L)', 'G2_Cl (mg/L)', 'ISS1 (mg/L)', 'PO4 (mg/L)', 'NH4 (mg/L)', 'NO3 (mg/L)', 'FE (mg/L)',
    'LDOM (mg/L)', 'RDOM (mg/L)', 'LPOM (mg/L)', 'RPOM (mg/L)', 'BG (mg/L)', 'DIAT (mg/L)', 'OTH (mg/L)', 'DO (mg/L)']
w2.multiplot(cwo_37, figsize=(15, 30), ylabels=ylabels)
w2.multiplot(cwo_75, figsize=(15, 30), ylabels=ylabels)

DWO (derived constituent outflows), Berlin Model, 2006

• SEG 37: dwo_37_wdo.csv
• SEG 75: dwo_75_wdo.csv

# Constituents list: Include an extra "constituent" to the extra column that will be added due to a trailing comma.
constituents = ['TKN', 'TP', 'CHLA', 'TSS']
dwo_37 = w2.read('BerlinMilton2006/dwo_37_wdo.csv', 2006, constituents)

# The following file has times only, no data
# dwo_75 = w2.TimeSeries('BerlinMilton2006/dwo_75_wdo.csv', 2006, constituents)

ylabels = ['TKN (mg/L)', 'TP (mg/L)', 'CHLA (mg/L)', 'TSS (mg/L)']
w2.multiplot(dwo_37, figsize=(15, 15), ylabels=ylabels)
display(dwo_37)

	TKN	TP	CHLA	TSS
Date
2006-01-01 00:00:00	0.0000	0.0000	0.0000	0.0000
2006-01-01 01:00:00	1.1375	0.0636	4.4195	7.2167
2006-01-01 02:00:00	1.1373	0.0636	4.4188	7.2096
2006-01-01 03:00:00	1.1370	0.0636	4.4178	7.1997
2006-01-01 04:00:00	1.1369	0.0636	4.4177	7.1983
...	...	...	...	...
2006-06-25 17:00:00	0.5503	0.0479	2.0518	0.7734
2006-06-25 18:00:00	0.5504	0.0479	2.0502	0.7729
2006-06-25 19:00:00	0.5504	0.0479	2.0493	0.7725
2006-06-25 20:00:00	0.5501	0.0479	2.0546	0.7726
2006-06-25 21:00:00	0.5495	0.0479	2.0687	0.7733

4222 rows × 4 columns